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The laser-excitation of Rydberg atoms in ultracold gases is often described assuming that the 
atomic motion is frozen during the excitation time. We show that this frozen gas approximation can 
break down for atoms that are held in optical lattices or microtraps. In particular, we show that the 
excitation dynamics is in general strongly affected by mechanical forces among the Rydberg atoms 
as well as the spread of the atomic wavepacket in the confining potential. This causes decoherence 
in the excitation dynamics - resulting in a dissipative blockade effect - that renders the Rydberg 
excitation inefficient even in the anti-blockade regime. For a strongly off-resonant laser excitation 
- usually considered in the context of Rydberg dressing - these motional effects compromise the 
applicability of the Born-Oppenheimer approximation. In particular, our results indicate that they 
can lead to decoherence also in the dressing regime. 

PACS numbers: 37.10.Jk, 32.80.Ee, 34.20.Cf 



Ultracold laser-driven Rydberg gases are a versatile 
platform to study the coherent quantum dynamics in 
strongly interacting many-body systems. One reason 
that makes these systems so appealing is the fact that the 
thermal energy is so low that the atoms can be considered 
to be frozen in place on the timescale of laser-excitation 
[Ul^. This absence of thermal atomic motion entails that 
the dynamics of ultracold Rydberg gases is entirely deter- 
mined by the competition of the coherent laser-excitation 
process and the strong interaction between the highly ex- 
cited atoms. This interplay results in an intricate and 
highly correlated excitation dynamics - its most promi- 
nent manifestation being the dipole-blockade: due to the 
strong interactions between Rydberg states only a single 
atom can be laser-excited within a certain exclusion vol- 
ume, which in turn gives rise to an enhancement of the 
laser coupling to the emerging many-body state p IJ. 
Both, the dipole-blockade as well as the enhanced laser- 
coupling are at the heart of possible applications of ultra- 
cold Rydberg gases in quantum information science [5], 
the simulation of quantum spin models ^61112]. and the 
creation of highly nonlinear and nonlocal optical media 

On a sufficiently long timescale dispersive forces among 
Rydberg atoms do eventually lead to atomic motion. 
However, this timescale is typically about one or two or- 
ders of magnitude larger than the excitation pulse dura- 
tion. Due to this mismatch an adiabatic approach is often 
used where one considers the laser excitation in the frozen 
gas limit and treats the atomic motion separately within 
a molecular dynamics framework [18^-^20]. A similar route 
is also followed in the theoretical description of Rydberg 
dressing protocols |2TH27i . where a strongly off-resonant 
laser coupling is used to admix a very small fraction of the 
Rydberg wavefunction to the atomic ground state. Here 
one conventionally applies a Born-Oppenheimer approx- 
imation, which consists of first determining effective po- 



tentials emerging from the off-resonant laser coupling for 
fixed atomic position and then solving for the motional 
dynamics of the Rydberg dressed system. 

In this work we show that for atoms which are trapped 
in optical lattices or microtraps the interplay between 
electronic and motional dynamics can already be impor- 
tant on typical timescales of resonant Rydberg excitation. 
To this end, we will investigate the quantum dynamics 
of a one dimensional model of the laser-excitation of two 
atoms, each trapped in a separate potential well. We 
will show that the excitation dynamics changes its char- 
acter from fully coherent to dissipative when the distance 
between the atoms is decreased. The dissipative dynam- 
ics emerges when the mechanical force between the Ry- 
dberg atoms is so strong that the adiabatic approxima- 
tion breaks down. This is accompanied with a signifi- 
cant slowing down of the excitation timescale leading to 
a dissipative blockade effect, where the laser excitation 
of a Rydberg atom pair gets inefficient even in the anti- 
blockade configuration. Moreover, our findings indicate 
this dissipation can also affect the dynamics of ground 
state atoms which are weakly dressed with a Rydberg 
state. In particular we will argue that the concept of an 
effective potential between the Rydberg dressed atoms 
can become insufficient to describe the dynamics of the 
system. Our results have implications on current and 
projected Rydberg experiments in optical lattices [281129] 
and microtrap setups [30-32 . 

In order to develop an understanding of the effect of 
atomic motion on the excitation process, we consider a 
model describing the laser excitation of a pair of atoms 
in one spatial dimension (cf. Fig. [T]). The electronic 
structure of each atom is modeled by two levels: the 
electronic ground state \g) and an highly excited state 
|e). These states are coupled by a laser field which is 
detuned by 6 from the atomic transition and we param- 
eterize the coupling strength by the Rabi frequency Qq. 
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FIG. 1. (a) Anti-blockade configuration. Two atoms in their 
electronic ground state are prepared at a distance ro and sub- 
sequently resonantly photo-excited to a Rydberg pair state 
with Rabi frequency Q. Each atom is initially assumed to be 
in the motional ground state of a harmonic potential (inset). 
The laser couples the corresponding wavefunction for the rel- 
ative coordinate (blue) to a continuum of unbound Rydberg 
pair states (red). The overlap of these wavefunctions deter- 
mines the energy-dependent laser coupling W{E), which is 
shown for two different values of the initial separation ro . For 
our analytical calculations we linearize the interaction poten- 
tial around r = ro (dashed line). The numerically obtained 
time evolution of the survival probability p{t) of the initial 
state is shown for the cases of large (b) and small (c) initial 
separation ro of the atoms taking into account the full po- 
tential V{r) = Ce/r^ for the Rydberg 65 S state of rubidium 
atoms with Ce ~ 27r x 370 GHz-/im^. The trapping poten- 
tial is cjt = 27r X 13.4 kHz yielding Ar ^ 133 nm. The Rabi 
frequency is set to be Q/oJt = 50. The initial separation is 
ro/Ar = 72 in (b) and ro/Ar = 30 in (c). The time interval 
shown corresponds to 0.96 /is. 



We assume that two ground state atoms are initially pre- 
pared at a distance ro, each in the lowest motional state 
of separate harmonic traps [see inset of Fig.jlja)]. This 
reflects the situation which is encountered in the case of 
atoms in a deep optical lattice or cooled to the motional 
ground state of microtraps [33l ^ . In this case the spa- 
tial wavefunction of each atom is a Gaussian of width 
a. For ro ^ cr the wavefunction describing the state 
of the relative motion of the atom pair is also a Gaus- 
sian x{r;ro) = (27r)"^/'^Ar"^/^exp [-(r - ro)^/(4Ar^)] 
of width Ar = V2a. 

Anti-blockade configuration - Our first goal is to study 
the dynamics of the direct laser-excitation of a pair of 
Rydberg atoms as depicted in Fig. [TJa). To achieve this 
the laser frequency needs to be chosen such that the anti- 
blockade condition A(ro) = 2S-\-V{ro) = is met [35l|36], 
where V(r) denotes the interaction potential of the Ry- 
dberg atoms. If |l^(ro)| is much larger than the Rabi 
frequency of the atomic transition, the anti-blockade 
condition also entails |^| ^ (^o- In this regime the pair 
states with only one Rydberg excitation can be adiabati- 



cally eliminated and the Rabi frequency of the transition 
between the pair states \gg) and \ee) having zero and two 
Rydberg excitations, respectively, is ft = ^^o/l^l- These 
considerations, however, do not account for the spatial 
uncertainty of the initial pair state caused by the zero- 
point motion of the ground state atoms in their potential 
wells. Since the anti-blockade condition is strictly valid 
only for r = ro, an uncertainty Ar in the atomic separa- 
tion will affect the laser excitation of the atom pair. In 
order to assess the consequences of this, we first investi- 
gate the dynamics of our model by solving it numerically. 

The Hamiltonian of our system can be divided into a 
part describing the atomic motion and an electronic part 
accounting for the laser-excitation, H = iJmot + ^ei, with 

{h = i) 

Hmot = -^V^ + U{r) \gg) {gg\ + V{r) \ee) (ee| (1) 
ffei = -A(r) |ee) {ee| + | {\gg) {ee| + h.c.) . (2) 

Here, m denotes the reduced mass of the atom pair and 
U{r) = muj^{r — ro)^/2 the harmonic trapping potential 
of the ground state atom pair with trapping frequency 
Ut. For simplicity we assume that the Rydberg atoms 
do not feel any external potential. In the following, we 
will focus on timescales that are much shorter than the 
radiative lifetime of the Rydberg atoms. 

Panels (b) and (c) of Fig. [l] show the numerically ob- 
tained [37 survival probability p(t) = | (G|e"^^^|G) |^ of 
the initial state \G) = |^^) 0x(^; ^o) as a function of time 
for large and small initial separation ro of the atoms, re- 
spectively. Evidently, the dynamics in these two regimes 
has a strikingly different character. While for a large 
separation (b) the initial state gets depopulated quickly 
and coherently, panel (c) displays an exponential decay 
of p{t). The major difference of these two regimes is 
the variation of the interaction potential V{r) over the 
spread Ar of the relative wavefunction. This implies 
that the laser actually does not just couple two discrete 
electronic states {\gg) ^ |ee)) but instead the discrete 
state \G) and a continuum of states \E) = \ee) (g) (j){r^E) 
of energy E within an energy window |F|Ar. Here 
F = — drV{r)\^^^^ is the force between the Rydberg 
atoms at distance ro- For a weak force the laser coupling 
is essentially constant over Ar and only continuum states 
within a very small energy window are involved in the ex- 
citation dynamics. In the case of a strong force, however, 
V{r) varies significantly over Ar so that the laser cou- 
pling is smeared out over a large energy interval. 

In order to obtain an analytical understanding of 
the excitation dynamics we employ the framework of 
Fano theory [38l We will assume in the fol- 

lowing that the oscillator ground state is energeti- 
cally well isolated from higher oscillator levels [40]. 
In this regime the Hamiltonian of our model takes 
on the form H = A{ro)\G){G\ + JdEE\E){E\ + 
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FIG. 2. (a) Dependence of the spectral function |a(c(;)|^ on the 
initial separation ro of the atom wavepackets for A(ro) = 
(see Fig. [T]) and Q/cJt = 1. The dashed lines indicate the 
positions of the eigenstates of the model when motional ef- 
fects are not taken into account. The distance Rc, where the 
spectrum changes from a single to a double peak structure 
is indicated by the vertical dotted line. Panels (b) and (c): 
Decay rate in the strong force regime as function of the ini- 
tial atom separation for Q/cJt = 1 (b) and Q/cjt = 50 (c). 
Comparison of the numerical solution of the model (symbols) 
and the analytical prediction obtained with Fano theory (solid 
line), (d) Decay rate in the strong force regime as function of 
A(ro) for two different values of the initial atomic separation 
and Q/uJt = 1. The symbols and solid lines are the numeri- 
cal solution and the predictions of Fano theory, respectively. 
Numerical parameters that are not explicitly indicated in the 
panels are those of Fig. [l] 



/ dEW{E) {\E){G\ +h.c.). The energy-dependent cou- 
pling between the discrete state and the continuum, de- 
noted by W{E), is proportional to the spatial over- 
lap between x(^5^o) with the continuum state (j){r^E) 
and given by W{E) = (Sl/2) J^drx{r;ro)^{r,E). The 
eigenstates of H can be expressed as = a{uj) \G) + 
J^^dE l3{uj,E)\E) with coefficients a{uj) and (3{uj,E) 
that have to be determined self-consistently [39] . Choos- 
ing the initial state to be \G) the probability to remain 

in it at time t then is p{t) = J^^du |a(co')|^ e~*^ 

The dynamics of the system is therefore encoded in the 
Fourier transform of the spectral function, |Q^(a;)|^ = 
W^{uj)/{[uj-A{ro)-e{uj)]^^7r^W'^{uj)}, with level- 
shift function e(cj) = V J^^ dEW'^{E)/{uj-E), where V 
denotes the principal value integral. The spectral func- 
tion is normalized such that do; |q^(cj)|^ = 1. 

To determine the spectral function analytically, we 
make two approximations: (i) We linearize the in- 
teraction potential around ro, i.e., we approximate 
V{r) ~ V{ro) + F{r — ro). The continuum states 



0(r, E) are then the eigenstates of a particle in a linear 
potential, ^{r,E) = M'^/^^ {- [{r - ro) ^ E/F] /lo}, 
where ^{x) denotes the Airy function [41 . The con- 
stant A/" = 7r|F|^/^/(2m)^ is chosen such that these 
eigenfunctions are normalized to a J-function in en- 

1 /3 

ergy and Iq = [l/{2m\F\)] ' . (ii) We consider the 
regime where Ar/lo 1. Within these approxima- 
tions the coupling acquires a simple form, W{E) = 

2^1/4 A£;i/2 exp (-2^)^ with energy width = 

^/2Ar\F\. In this Gaussian coupling regime the level 
shift function can also be expressed analytically, e{uj) — 

^Ke^ {'Ke) ' where V{x) = dy is the Dawson 

function. 

The spectral function |a(co')|^ in the anti-blockade 
configuration [A(ro) = 0] is shown in Fig. |2|a) as a 
function of the initial distance of the atomic pair for 
Q/ujt = 1. It displays a transition from a single peak 
to a double peak structure at a critical distance Rc = 
1.363[y(Ar)/r^)]^/^Ar assuming a van der Waals inter- 
action, V{r) = C^jr^ ^ between the Rydberg atoms. In 
this transition region the spectral function is broad, while 
for small/large ro, i.e., strong/weak force, its peaks are 
very narrow. For ro/Ar ^ 1 the peaks are located at = 
±1^/2. This is the result one would obtain for the anti- 
blockade irrespective of the atomic separation if one com- 
pletely neglected the variation of the interaction potential 
over the spatial width of the initial wavefunction. In this 
case, which corresponds to panel (b) in Fig. [l] the spec- 
tral function is |a(a;)|^ = [(5 (cj + 1^/2) + 5 (cj - 1^/2)] /2. 
It has two (5-peaks at the eigenenergies of the electronic 
Hamiltonian ([2| which are the dressed energies of the 
fully coherent system [depicted also as dashed lines in 
Fig. |2|a)]. The change of the spectral function with the 
atom separation for finite Ar therefore clearly shows that 
motional effects can have a significant effect on the exci- 
tation dynamics when the force between Rydberg atoms 
is strong. 

To study this further, we focus on the regime, where 
ro <C i^c- Here, the coupling WiE) is spread over a very 
large energy interval A^, such that it can be considered 
constant as a first approximation. The spectral function 
in this regime acquires a Breit-Wigner form. 



r/2 



7r(r/2)2 + {c.-A(ro)-e[A(ro)]P 



(3) 



i.e., the survival probability of the initial state decays 
exponentially, pit) = exp(— Ft), as shown in Fig. [TJc) 
with rate 



2 AE 



exp 



A^(ro) 



(4) 



Let us compare these predictions with results obtained 
by numerically solving our model and first focus on the 
regime, where the assumption of an energetically well iso- 



4 



lated initial state is well met. Fig.[2[b) shows a compari- 
son of r with the decay constant extracted from fully nu- 
merical simulations as a function of the initial atomic dis- 
tance for A(ro) = 0. The numerical data were obtained 
by analyzing the time evolution of the survival proba- 
bility of the initial state (To) and by following the dy- 
namics of the total population in the harmonic potential 
(Ftot)- The analytical predictions are in excellent agree- 
ment with the numerical data. In addition, also for finite 
detuning the agreement between the numerical simula- 
tions and the analytical prediction of Eq. Q is very good. 
This can be seen in Fig.[2|d) showing the rate as function 
of A(ro) for two different initial separations. The slight 
asymmetry of the numerical data about A(ro) =0 stems 
from the fact that the force between the Rydberg atoms 
is not constant over the initial wavepacket as assumed 
in our analytical calculations. Interestingly, as demon- 
strated in Fig.|2|c), our analytical analysis also gives rea- 
sonable agreement when the assumption of an energeti- 
cally isolated initial state does not hold. At small spatial 
separations the agreement between numerical and ana- 
lytical results is remarkably good. For larger ro/Ar the 
full solution shows that the oscillator ground state gets 
coupled to higher levels via the continuum. However, ir- 
respective of this, the excitation dynamics dramatically 
slows down with increasing mechanical force between the 
Rydberg atoms and the excitation rate is inversely pro- 
portional to it. Thus, although on resonance, the Ryd- 
berg excitation becomes more and more inefficient with 
increasing mechanical force. This effect can be viewed as 
a dissipative excitation blockade induced by decoherence 
due to atomic motion. 

Dressing regime - Let us finally analyze the situation, 
in which the excitation laser is far-detuned from the sin- 
gle atom transition as well as from the anti-blockade con- 
dition. This is the regime of Rydberg dressing, where 
only a very small fraction of the Rydberg wavefunction 
is admixed to the ground state atoms [HI |22l . The 
main effect of this small admixture is that a pair of Ry- 
dberg dressed atoms exhibits a distance-dependent en- 
ergy shift. In the Born-Oppenheimer approximation, the 
spectral function consists of two (5-peaks at the eigenen- 
ergies of the r-dependent electronic Hamiltonian - as in 
the anti-blockade configuration. The lowest eigenenergy 
corresponds to the dressing potential and is given by 
uj{r) = A(r) - v'A2(r) + (^2 /2. The black dashed line 

in Fig. [3]^ a) depicts the interaction energy of two dressed 
atoms in the motional state x(^5^o) which is given by 
the convolution uJoi'^o) = / dr \x{r;ro)\'^uj{r). For a van 
der Waals interaction this energy exhibits a character- 
istic soft-core with radius = {C6/2\S\)^^^ and height 
Vq = ^}q/8\S\^. Moreover, the probability to be in the Ry- 
dberg pair state is given by, Pee(^o) = / dr |cee(r, ro)P, 
where Cee(r, ro) is the coefficient of the Rydberg pair state 
in the eigenstate \oj{r)) = Cgg{r,ro)\gg) + Cee(r, ro)|ee) 




FIG. 3. (a) Distance dependence of the energy (black dashed 
line) of a Rydberg dressed atom pair (jJD{ro) in the motional 
state xi'^^'^o)- (b) The spectral function of the sys- 

tem showing a significant spectral broadening as the soft-core 
radius is approached. The peak position of |a(cj)|^ is shown 
as solid red line in panel (a) up to the point where the width 
of the spectral peak starts to exceed its height, (c) Distance 
dependence of the probability to be in the Rydberg pair state 
for t ^ 00. The data set obtained from the Fano theory 
is given by the circles while the solid line shows the results 
from the Born-Oppenheimer approximation. The numerical 
parameters used are ^ = 27r x 20 MHz, Qo = 27r x 1 MHz and 
rs/Ar ^ 34.5. 



corresponding to the eigenvalue Lj{r). This probability is 
shown as a function of ro as solid line in Fig. [sj^c). 

This picture is changed when the effects of atomic mo- 
tion are included in the description. Fig. [sj^b) shows the 
spectral function |a(cj)p obtained from Fano theory. For 
ro ^ 00 one recovers the S-peak at the eigenenergy of 
the dressed system. However, as the soft-core radius r^ is 
approached this peak broadens and the spectral weight 
is spread over an energy interval largely exceeding Vq. 
Eventually the broadening is so large that the spectral 
function is almost ffat and no pronounced structure is vis- 
ible in |a(cj) p. Moreover, as shown by the solid red line in 
Fig. [3]^a) the peak position is shifted as compared to the 
adiabatic energy of the dressed state. We show the data 
up to the point where the peak height becomes smaller 
than its width. The large broadening of the spectral func- 
tion entails that the energy of the Rydberg dressed state 
gets less and less well defined with decreasing atomic sep- 
aration. This suggests the emergence of a dissipative dy- 
namics for sufficiently small interatomic distances. Com- 
paring the excitation probability to be in the Rydberg 
pair state in the limit t ^ 00, Peei^o) = f dcj lafcj)!^ 
(the integration limits (c^min/max) are chosen such that 
they contain the entire peak), we find that this prob- 
ability actually coincides with the results of the Born- 
Oppenheimer approximation [see Fig.|3|c)]. 

Neglecting the effect of atomic motion, one thus ob- 
tains the correct probability to be in the electronic states 
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\gg) or |ee), however, without further information about 
the character of the excitation dynamics. This ihustrates 
that the inclusion of motional effects in the description 
of the excitation process can be important even in the 
far off-resonant regime of Rydberg dressing. 
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